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Abstract 



O . 

(yQ , We discuss the multiplicity distribution for highest accessible energies of 

G\ , 

i pp- and pp- interactions from the point of view of the multiparton collisions. 

O-t'. 

' i The inelastic cross sections for the single, ai, and multiple (double and, pre- 

(D , 

r~| , sumably, triple, (T2+3) parton collisions are extracted from the analysis of the 

, experimental data on the multiplicity distribution up to the Tevatron ener- 

^ , gies. It follows that ai becomes energy independent while (T2+3 increases with 

^/s for ^/s > 200 GeV. The observed growth of {p±) with multiplicity is at- 
tributed to the increasing role of multiparton collisions for the high energy 

pp{pp)- inelastic interactions. 
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I. INTRODUCTION. KNO SCALING AND ITS VIOLATION 



The experimental observations of the violation of KNO scaling |T| at the high energy 
(pp)- and (pp) -collisions and of the correlations between the average transverse momen- 
tum {p±) and the multiplicity of the secondaries (a higher {p±) for high multiplicity 
events) [^] indicate that there are at least two (complementary) mechanisms of the high 
energy multiparticle production revealing the quark-gluon structure of hadrons and their 
interactions. 

That KNO scaling should be violated at very high energies was realized long ago ^ , in the 
same year when KNO scaling was introduced. This deviation from a KNO type distribution 
of secondaries was attributed to the possibility of the splitting of each of the colliding 
hadrons into several constituents (valence quarks, partons) pair wise interacting with their 
counterparts from oppositely moving hadrons. The above picture results in the production 
of several showers and, in the frame of the Reggeism, is attributed to the contribution of 
Regge cuts to the elastic scattering amplitude, in addition to the Regge pole (Pomeron) 
corresponding to the production of the single shower. 

Inclusion of such contributions changes the structure of the distribution (aM is the 
cross section of the production of secondary hadrons) at large ^ = In ^ (sq = 1 GeV^) 
leading to the appearance of the additional peaks in the distribution with larger (A^„) (here 
n denotes the number of pairs of the simultaneously colliding partons involved in the in- 
teraction from the different hadrons, or the number of resulting showers; n = 1, a single 
interaction, can be attributed to the Regge pole; n = 2, a double collision of two pairs of 
partons from different hadrons, can be attributed to the exchange of two Pomerons (Regge 
cut), etc.). 
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II. REGGE PICTURE AND MULTIPLE COLLISIONS OF PARTONS 

In this picture, (Nn) = n{Ni) where (A^i) is the mean multiphcity corresponding to 
the production of one shower in the single partonic colhsion. We have, from the data 
where the one Pomeron exchange describes experiments at lower energies (A^^i) = a + 61n^, 
with a = —3.879; b = 2.099 (see Appendix A for details). This qualitative picture 
corresponds to the idealized scenario where only "non-enhanced" Regge-type diagrams are 
included ( "quasi-eikonal" approximation). Inclusion of the so-called "enhanced" diagrams 
where the possibility of the branching of ladders-showers is allowed, thus giving the emission 
of the smaller average number of emitted particles than the corresponding non-enhanced cut, 
leads to the smearing out of the distribution between the peaks reducing the effects of KNO 
violations. Since the distributions of particles in the different ladders are independent, one 
expects that AA^„ ~ y/nNi giving the broadening of the distribution 

We emphasize that the multipomeron exchanges are especially important in the modern 
treatment of the Pomeron as a pole with intercept a{0) higher than unity (a(0) = 1 -|- 
A; for "soft" Pomeron A ^ 0.08 as experimental data on high energy "soft" collisions 
show), where the relative role of these exchanges is increased with energy (e.g., Pomeron 
contribution to the o"tot is ~ e^^, whereas the n Pomeron exchange is ~ e"^^). Concerning 
the effective number of exchanges in "quasi-eikonal" approximation, one has ries ~ e^"^ < 
2-5 (^^n = ^^ij^)^ Hi "\A ^ 2 TeV (n = 2). Thus, from this point of view, only double 
parton interactions (two effective Pomeron exchanges) in addition to the single collision are 
expected to be effective at the energies up to the Tevatron energy. At LHC energy {^/s = 14 
TeV) this parameter ~ 3.3 {n = 3), thus indicating the possibility of the appearance of 
the third maximum in distribution with {N^) ^ 3(A'i) (see Fig. 1 and Appendix B). 



^ Since at the highest achievable energies {^/s > 500 GeV) KNO scaling (in the whole rapidity 
interval) is violated strongly, this indicates that at these energies the "enhanced" graphs, probably, 
are not important. It is not excluded, that at higher energies KNO scaling will be restored. 
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III. GLUON STRING MODEL AND POMERONS 



On the QCD level, the ultra-high energy interaction between hadrons phenomenologically 
described by Pomeron exchange is treated as a result of a (few) gluons exchange between the 
constituents of hadrons as they pass close to one another [0. After a color exchange between 
partons (valence quarks, gluons, etc.) carried by the colliding hadrons, partons from one 
hadron are joined by two gluon strings with partons inside the other hadron. This pair 
of strings stretch when hadrons separate after the collision and break down transforming 
into the system of hadrons forming a Pomeron-like shower (ladder), giving the asymptotic 
density of hadrons. 

From this point, the above described n Pomeron exchanges are considered as the joining 
and stretching of n different pairs of strings between hadron constituents resulting, after the 
hadron separation, the production of n Pomeron-like showers with n-fold increased average 
multiphcity (A"™). 

This picture using the language of gluon strings can be related to the Pomeron phe- 
nomenology, at least in the quasi-eikonal approximation, and, as shown in leads to the 
same expressions for a variety of characteristics of the high energy strong interactions as the 
Regge diagram technique, giving also the possibility to determine some free parameters of 
last one. 

IV. TWO COMPONENT PICTURE OF THE HADRON COLLISIONS 

From the more phenomenological approach [^j, it is possible to argue that there are at 
least two components in the cross sections for the ultra-fast hadron collisions: soft component 
corresponding to the single collisions of partons from different hadrons which can (but not 
necessarily) be described by the exchange of the single Pomeron (or the pair of gluonic 
strings), and the so-called semi-hard (jet) component which is responsible for the higher 
of the secondaries than is usually attributed to the single "soft" Pomeron exchange. 
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At low energies (up to energy of ISR) only soft interactions with scaled multiplicity 
distribution of the KNO type are important. At higher energies (CERN SPS, Tevatron) 
the KNO-like shape of the multiplicity distribution is modified by the increasing role of 
the "hard" multiparton interactions ^ (four parton collisions of ||T0[) that leads to the 
increase of the fraction of events with higher multiplicity {N ^ (Ni)). This new channel 
is characterized, at the same time, by the higher {p±), in accordance with the experimental 



observation ||Tl| that the high multiplicity events indicate a jet-like structure. 

From the point of view of the interactions between constituents, one can argue that this 
second "semihard" component with (A^2) ~ 2(A'^i) is due to the double collisions when in the 
same inelastic interaction, two different parton pairs from colliding hadrons independently 
collide. Since the parton densities of the colliding hadrons overlap strongly at small impact 
parameters (higher p_i_) the probability of the double parton collision {n = 2) is enhanced, 
and, as a result, two pairs of soft partons (e.g., valence quarks) {x± = ^ <^ 1) stacked in 
the collision resulting in the (mini) jet structure of the secondaries with > (A"2) = 2{Ni). 

In 1^ it was shown that the above two-component model gives a successful fit for the total 
non-single diffractive inelastic cross section inclusive spectra and multiplicity distribution in 
full phase space for the energies y/s = 200 and 546 GeV. We note that in ^ the ad hoc 
assumption was made that the inelastic cross section for the soft component (single parton 
collisions) is energy independent. 

We will see below, from our analysis of the multiplicity distribution, that the high energy 
inelastic cross section corresponding to the single parton collisions indeed, practically does 
not depend on y/s, as it was assumed in 0. On the other hand, the one Pomeron exchange in 
its quasi-eikonal form gives an increase of inelastic cross section. We don't know the answer 
to that inconsistency, except that the role of the enhanced diagrams is underestimated, and 
they have to be included in the analysis of the "soft" cross section. 
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V. VERSUS MULTIPLICITY 



Now, if we again turn to the language of the Pomeron exchanges, it is possible to con- 
jecture that the above described double parton collisions can be attributed to the double 
effective Pomeron exchange which leads not only to the higher multiplicity (and to the devi- 
ation from KNO shape) but also to the increase of the average transverse momentum {p±) in 
the regime of ^ <^ 1. Indeed, the p_i_ dependence of n Pomeron exchange amplitude is given 
by the expression exp ^— ''^"^-^ ^ , where A„ = + a'(0)^„. Experimentally, for pp-coUisions, 
a'{0) = 0.25 GeV-2 g. That gives for the ratio 



This ratio is in good correspondence to the experimental data |T^. The further increase of 
the {p±) can be expected at the effective opening of the triple parton collision, which as we 
see from Fig. 1 should be clearly displayed at the LHC energy range. An experimentally 
established increase of {p±) with energy is connected, of course, with the gradual increase of 
the relative role of the double Pomeron exchange (or double parton collisions) with energy. 

There are recent results from the CDF collaboration |jl3[ for the process — independent 
and free from theoretical schemes cross section of double parton collision in the 3 jets data: 

+1.7 

o-dp,eff = (14.5 ± 1-7 _2 g) mb. 



The AFS experiment |jTj| (y^ = 63 GeV) gives for the same quantity ~ 5 mb, while UA2 
presents a lower hmit of ddp^es > 8.3 mb |jl5| . 



VI. ANALYSIS OF THE MULTIPLICITY DISTRIBUTION DATA 

Here we analyze the data on the multiplicity distributions obtained at the energy range 
^/s from 30 to 1800 GeV including ISR, UA5 and Tevatron (E735 experiment) data. We use 
as a basis of our analysis the fact that KNO scaling is well satisfied in the experimental data 
through the ISR energies. The deviation from the simple KNO scaling at higher energies is 
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due to another process which is incoherently superimposed on the KNO producing process. 
In Fig. 3, the differential cross sections are displayed for the UA5 and the E735 data. Note 
that the E735 data is likely to be less reliable at low multiplicities than the UA5 data whereas 
the E735 data is more statistically rehable at the higher multiphcities. 

In Fig. 3 the coUider data from different energies is superimposed on one plot. The data 
have been normalized to the same value of the variable x = at the Xj^ax, maximum 
value ttsd^- It turns out that the cross section at the maximum is essentially indepen- 
dent of energy. From the measurement of KNO curve at the ISR energies we know that 
-/Vmax = 0.8(A^i). Thus, we can estimate the value of (iVi), the average multiplicity for the 
processes with pure KNO scaling from the known value of A^max at which the cross section is 
a maximum for higher energies. We made a least square fits to the value of Amax and (A^i) 
over range ^/s from 30 GeV to 1800 GeV for the former and ^/s from 30 GeV to 200 GeV 
(where KNO scaling is well satisfied) for the latter. These two quantities are well fitted by 
the form a + 6 In i/s with the same coefficients over their regions of overlap. As we remarked 
above, the deviation from KNO scaling is due to another process which is incoherently su- 
perimposed on top of the KNO producing process. By subtracting the KNO distribution 
from the observed data we determine the shape of the competing process as shown in Fig. 
4. 

The shapes of the multiplicity distribution thus found are rather different from the KNO 
distrubition and is not a simple convolution of the KNO distribution. There must be a 
correlation between the multiplicities in the two collisions. The main characteristics of the 
derived distribution is that the most probable value of the distribution occurs at a; = 2 (or 
at twice the multiplicity corresponding to the initial low energy KNO distribution). The 
width of the distribution is close to \/2 times the width of that KNO shape. This is in quite 
good agreement with the picture we presented above, which is based on the adding of the 
double collision of partons of the colliding hadrons by the exchange of the two Pomerons (or 
two pairs of strings) to the single parton collision which proceeds by the exchange of one 
effective Pomeron. 
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Anyway, independently of the concrete realization of the interaction between ultrarela- 
tivistic parton pairs, we can interpret the population of the secondaries with maximum at 
2(A^i) as a result of two independent parton-parton collisions occuring in the same encounter. 
We denote the inelastic (non-single diffractive) cross section for this process by 02 whereas 
(Ti denoted the inelastic cross section due to the single parton-parton collision. This cti is 
characterized by the KNO scaling and, as it follows from our analysis, it is nearly energy 
independent for ^/s > 200 GeV and has a value of (34 ± 2) mb (Fig. 5). 

Some remarks on the experimental data of E735 and UA5 are necessary. In Figs. 3-5 we 
do not show error bars on all the points in the derived distributions. The statistical errors 
in E735 are relatively small. However, the systematic error in both experiments might be 
sizeable. In E735 as well as UA5 the multiplicities in a restricted range of rapidity are 
extended to the full range by computer simulation. There are corrections for secondary 
interactions in individual events. The derived distributions in Fig. 3 suffer in accuracy from 
the fact that one uses subtraction of two different distributions which each have uncertainties. 
When the experimental data and the KNO distribution are close in value then the error is 
greater. Fortunately, the KNO distribution is only the order of 10% of the experimental 
data at = 2(A^i) [x — 2) and consequently, the position of the peak in the derived 
distribution is not strongly affected by errors produced by taking the differences. The 
whole distribution might move up and down as a result of errors in either of the KNO or 
experimental distrubitions. 

The values of the cross sections represented by the curves shown in Figs. 3,4 are given 
in Fig. 5. 

It is seen there that, in contrast of cti, (T2 is increased with ^/s and equals 16 mb at 1800 
GeV. We already have mentioned that CDF collaboration gives the effective double-parton 
collision cross section equal (14.5 ± 1.1^^''^) mb. One has to reahze that in the sample of 
events with double parton coUisions CDF has removed events with possible triple collisions. 

Wc can conclude that the double (and possibly, triple) collisions account for a large 
fraction of the increase in the total pp-cross section and, definitely, are responsible for the 
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increase of the pp inelastic cross section. 

As the coUision energies are increased to LHC values, it seems likely that double and 
possibly triple collisions will constitute a larger fraction of the inelastic cross section as is 
seen in our effective analyses from the previous sections (see Fig. 1). 

Another indication of the increasing importance of double collisions comes from an ex- 
amination of the (p±) of the charged particles at high multiplicity events. As we mentioned 
above there is an increase of {p±) as multiplicity increases. We note that {p±) increases by 
close to a factor from low to high multiphcity. The picture with the effective double 
Pomeron exchange in the double-parton collision leads just to this \/2 extra factor. The 
increse of the fraction of the multi-collisions with energy will, in general, lead to this well 
known (p^) -multiplicity relation. We expect that at LHC energies {p±) will undergo the 
additional increase. 

VII. CONCLUSIONS 

In this paper we have attempted to isolate and study the double-parton collision mech- 
anism from the analysis of high energy multiplicity data. The main result of our analysis 
of these data is that the non-single diffractive inelastic cross section consists of two compo- 
nents. The first component corresponding to the single parton collisions is characterized by 
being practically independent of ^/s cross section for ^/s > 200 GeV, whereas the second 
part of the cr^^^ increases significally with energy and achieves the value 16 mb at ^/s = 1.8 
TeV. That part was attributed here to the inclusion in the collision process of the double 
(and, maybe, triple) parton collisions. Thus, the increase of the aj^^^ at high ^/s is almost 
entirely due to the multi-parton collisions. 

On the other hand, we know from the experiment at lower energies where the single 
parton collisions indeed dominate that <J^^^ is increasing with ^/s. This indicates that at 
higher energies ^/s > 200 GeV the inelastic cross section due to the single coUision goes to 
the saturation, whereas the double coUisions give the increasing part of the cr^^^. One can 



9 



go further and conjecture that the same saturation will occur for the part of (7i„ connected 
with the double parton collisions at much higher energies whereas the (Tj„ due to the triple 
collisions will still increase with ^/s, etc., until asymptotically the total inelastic cross section 

(without the diffractive part) will set a constant value Gin — o"!^^*'* + (72^'^*^ + Of course, 

we cannot say anything about the behavior of Uei and Udif at this hypothetical limit. 

We need to emphasize that in this paper we have often used the notion of the "soft" 
Pomeron exchange in the process of parton collisions. Of course, this notion should be 
understood as an effective tool of the description of the ultrarelativistic parton collision 
and have played only the illustrative role and the trend in our analysis. However, the 
picture of the multi-parton collisions described here and the corresponding values for ai and 
(T2+3 extracted from our above analysis of the multiplicity are independent of the concrete 
mechanism of the parton-parton interaction. 
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FIGURES 

FIG. 1. a) Topological cross sections gn in the quasi-eikonal approximation with exchanges 
of three effective "soft" Pomerons for = 546, 900, 1800 and 14.10'^ GeV. b) Topological cross 
sections resulting from double and triple parton collisions for ^/s = 546, 900, 1800 and 14.10^ GeV. 

FIG. 2. Same in Fig. 1 for ^ = 100, 200 and 300 GeV. 

FIG. 3. A comparison of multiplicity distributions at different values y/s. The distributions are 
normalized at the maximum value of ^ where x = The solid curve is the KNO distribution 
from ISR data. 

FIG. 4. Multiplicity distributions obtained by taking the difference between the pp collider 
data and the KNO distribution. 

FIG. 5. Cross sections for the single (cri) and multiparton (double and triple) (0-2) parton 
collisions clS cl function of a/s- 
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APPENDIX A 



According to our ideology, at high energies {^/s ^ 1 GeV) the single parton collisions 
which may be described by the effective single Pomeron exchange and which are character- 
ized by the scaled multiplicity distribution, must be applied only up to the energy range 
where KNO scaling takes place. That means that the average multiplicity of the hadrons (in 
pp- and pp- interactions) for the single collisions (A^^i) must be obtained only from this energy 
range. The usually used total range of the existing energies up to highest ones (Tevatron), 
being important, of course, however, from the developed point of view, disguises the mech- 
anisms of the high energy collisions. Thus, we describe, in the spirit of the Regge picture, 
the average multiplicity (A^^i) as a + 61n ^ and find the coefficients a and b from the data of 
the high energy pp-collisions in the range of ^/s (11.5 ^ 62, 6) GeV [|161Q As a result of the 
fit, we have for the coefficients a and b: a = —3.874, b = 2.099 with = 7.085/7. 

We use the corresponding (A^i) as input in calculation of the multiplicity distributions 
which takes into account the multiparton {n = 2,3) collisions (see Appendix B). 

APPENDIX B 

We briefly describe here the results for the cross section of the production of particle in 
the quasi-eikonal approximation of the Regge approach (or, equivalently, as we emphasized in 
the text, in the model of the gluonic strings [^||). These results are represented graphically 
in Fig. 1. 

We are interested in the inelastic cross-section ain{N,s) not including the diffraction 
processes. The ignorance of the "enhanced" graphs allows us to consider that the long 
distant correlations among particles belonging to one shower-ladder are absent ^ . Similarly, 

^When our paper was written we became aware of the recent similar treatment of the average 
multiplicity (Ni) as a linear logarithmic function for "soft" component 
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we suppose that the partons from different showers are uncorrelated also. This permits us 
to assume that the distributions Pn{N) of the N produced particle (in full phase space) in 
the n shower events is the Poisson type:0 

Pn{N) = ^ e-<^") (1) 

{Nn) = n{Ni) =n(a + b\n—) . (2) 

We can write for ain{N, s) 

a,n{N, s) = a,P,{N) + a2{P2{N) + a,Ps{N) + ..., (3) 

where ct„(^„) is the cross section for the production of n ladders and ^„ = In (^^)- For 
we have [|] 



with 



(Tp = 87r7 — , 



The parameters in these expressions are fixed by the fits of the experimental data on fXtot 
and for pp- and pp-coUisions : 

7p = 3.64 (GeV)~^ = 3.56 (GeV)"^ C = 1.5, 
A = ap{0) - 1 = 0.08, a'p = 0.25(GeV)"2(c = 1). 



^Use of the concrete form of (1) is not crucial for obtaining the relation (5). It has a place for 
normalized Pn{N) with {Nn) = n{Ni) and large (A'^i). The relation (4), of course, does not depend 
on (1). Figs. 1-2 show only the trends. 
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We retain three terms in since the quadruple (or four-parton pair) coUision has small 
effect (see Fig. 1 which shows that even the triple collisions appear noticable only at the 
energe range of LHC). 

Summing (|^) over from 1 to oo with (Nn) ^ 1 we obtain for the total inelastic cross 
section 

CTinis) = CTi + 0-2 + as + 0((T4) (5) 

Fig. la) shows the A^-distributions of ain{N,s) at y/s = 0.55; 0.9; 1.8 and 14 TeV. At 
lower a/s the second peak of o"(A^, s) is not resolved. Fig. lb) shows the behavior of the 
parts of <j{N, s) which can be attributed to the double (and triple) collisions. 

In the table we present values of the inelastic "partial" cross sections (Xi, (T2 and us 
corresponding to the single, double and triple parton collisions in the Regge quasi-eikonal 
approach. 

We see that in this approach all [i = 1, 2, 3) increase with ^/s whereas our analysis 
in the text above indicates that the "single" collision contribution ai extracted from the 
experimental data on the multiplicity distribution practically is independent of ^/s. 

Furthermore, (Xi + (T2 + (T3 are systematically lower than the experimental cross section 
cr*°*(s) for the non-single diffraction events, whereas the corresponding theoretical values of 
this cross section resulting from the summation of all quasi-eikonal pomeron graphs 

where 

m - E 4^ - 1 Z <^ - ^"^T - -z (^c- ^' + '■'(^-^'' • 

with '-jE = 1-78 . . .. Euler constant and T{a, Z)- incomplete gamma function, are in excellent 
agreement with experimental data for cr*°*. 

Fig. 2a and 2b show for completeness the o"(A^, s) for lower ^/s (but higher than 62 GeV) 
{y/s = 100, 200, 300 GeV). On these figures at i/s = 200 GeV the shoulder, corresponding 
to the double collision is clearly seen. 
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TABLES 



V^, TeV 


(Ti, mb 


(72, mb 


(73, mb 


+ (72 + (73, mb 


0.55 


21.94 


9.57 


5.25 


36.76 


0.9 


22.72 


10.16 


5.71 


38.59 


1.8 


23.84 


10.72 


6.34 


40.90 


14 


27.19 


14.70 


8.43 


50.32 
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